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CN '. Abstract 

The dissociation of hydrogen molecules on the 7-U(100) surface is systematically studied with 
' the density functional theory method. Through potential energy surface calculations, we find that 

hydrogen molecules can dissociate without any barriers on the clean 7-U(100) surface. After careful 
electronic analysis, it is found that charge transfer between the hydrogen s and uranium d electronic 
^ ' states causes the dissociation, which is quite different from the dissociation of hydrogen molecules 

■ on other actinide metal surfaces. Considering that doping of 3(i transition metal atoms can stabilize 

O . the 7 phase of U, we also study the influences of Nb-doping on the hydrogen dissociation process. 

We find that the 2>d electronic states of Nb also take part in the hybridization with hydrogen s 
electronic states, which leads to the result that hydrogen molecules also dissociate without any 
energy barriers on the doped U surface. In addition, the free electronic energy lowers down more 
\ quickly for a hydrogen molecule approaching the doped U surface. 
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I. INTRODUCTION 



During the past decades, many theoretical and experimental studies have been carried 
out on the initial processes of surface hydrogenation or oxidation reactionsi*"— . The theo- 
retical researches include calculating the potential energy surfaces for H2 or O2 molecules 
on solid surfaces, and revealing the adsorption and dissociation mechanisms of thent^^— . In 
the experimental aspect, ultrahigh vacuum experiments are designed to detect the atomic 
structures after very few numbers of molecules reacting with solid surfaces^"". Actinides, 
as a group of radioactive, toxic and rare materials^, have their difficulties to be prepared 
for such experimental studies. On the other hand, they play important roles in advanced 
nuclear fuel cycles, and surface hydrogenation/oxidation are the main corrosion mechanisms 
that fail their storag Hence, theoretical studies are crucial for understanding the de- 

tailed surface corrosion mechanisms in the presence of environmental gases for these high-Z 
elements. Moreover, these studies may also help to address the environmental consequences 
of nuclear materials. 

Among the actinide elements, uranium (U) is the heaviest naturally occurring one, which 
occupies a central position in the early actinide series^i2i. Due to its important role in 
nuclear reactors, U is quite familiar to people^^. In the atmosphere, U and its alloys are 
ready to be oxidized to form U oxide layers, which subsequently break down through a 
hydrogenation process. The hydriding reaction proceeds by surface nucleation and growth 
of hydride nuclei which spread over the sample surfaces^^"— . It has been revealed that the 
hydride nuclei form beneath the oxide layers, and most of them are capable of penetrating 
through the oxide layers above them to grow over time until they consume the sample 
surface2Ii2^. Based on these knowledges, one can see that the electronic interaction with 
hydrogen molecules is critically important for the corrosion of U surfaces. Pure uranium 
crystallizes into several structures, the orthorhombic a phase with four molecules per unit 
cell at ambient conditions, followed by the body-centered tetragonal /3 (bet) phase at 940 
K and then the body-centered cubic 7 (bcc) phase at 1050 K at ambient pressure^i^i22i"— . 
Moreover, the high temperature 7 phase can be studied at normal temperatures by the 
addition of certain metals like molybdenum and niobium, which stabilizes the 7 phase at 
room temperature and below^i^ In metallic U, the three 5/ electrons of U hybridize with 
the 6(i and 7s electrons, and show itinerant behavior si^iS^i^^. Therefore, density functional 
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theories without any modifications concerning locahzed electronic states are appropriate 
within a large extent to describe metallic U. 

Recently, many theoretical and computational studies emerge on the surface chemical 
properties of uranium and other actinide metal surfaces. Using the density functional semi- 
core pseudopotential method, the chemisorption of carbon monoxide^ and oxygen gases^i on 
the 7-U surfaces have been investigated focusing on the geometric, magnetic and electronic 
properties of the system. The adsorption of carbon monoxide on a-U surfaces has also been 
studied using a plane-wave ultrasoft pseudopotential^^. For the interaction between hydro- 
gen and U, present theoretical studies are all carried out on the a-U surfaces^i^*^. Different 
from the hydrogen dissociation on transition-metal surfaces, where electronic hybridizations 
between hydrogen Is electrons and surface d electrons induce the dissociation^^, the dis- 
sociation of hydrogen molecules on the a-U surface is because of 5/-ls hybridizations^^, 
which is similar to the hydrogen dissociation on 6-Pu surfaces^I. Here in this paper, by 
using first-principles calculations, we systematically study the adsorption and dissociation 
processes of hydrogen molecules on the 7-U(100) surface, and reveal that the d electrons 
of U, instead of its / electrons, take part in the electronic hybridizations with hydrogen Is 
electrons. This result clearly presents different surface chemical properties of a- and 7-U. 
Considering that the 7 phase of U is available at normal temperatures only after doping 
with 3d transition metals, here we also investigate the influences of surface Nb doping on 
the dissociation properties of hydrogen molecules. 



II. CALCULATION METHOD 



Our calculations are performed within density functional theory using the spin- 
polarized edition of Vienna ab initio simulation package (VASP)^^. The Perdew-Burke- 
Ernzerhof (PBE)^'^ generalized gradient approximation and the projector-augmented wave 
potential^ are employed to describe the exchange-correlation energy and the electron-ion 
interaction, respectively. The cutoff energy for the plane-wave expansion is set to 520 eV. 
The U(IOO) surface is modeled by a slab composing of five atomic layers and a vacuum 
region of larger than 15 A . The bottom two layers are fixed, and the other U layers are free 
to move during the geometry optimizations of the uranium surfaces. A 2x2 supercell, in 
which each monolayer contains four U atoms is adopted in the study of the H2 adsorption. 
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Our test calculations have shown that the 2x2 supercell is sufficiently large to avoid the 
interaction between adjacent H2 molecules. Integration over the Brillouin zone is done using 
the Monkhorst-Pack scheme^ with 5x5x1 grid points. And a Fermi broadening^ of 0.2 
eV is chosen to smear the occupation of the bands around the Fermi energy Ep by a finite- 
T Fermi function and extrapolating to T=0 K. The H2 is placed on one side of the slab, 
namely, on the top surface. The calculation of the potential-energy surface is interpolated 
to 350 points with different bond length (dn-H) and height of the mass center [h-n^) of H2 
at each surface site. The calculated lattice constant of 7-U and bond length of a free H2 
molecule are 3.43 and 0.75 A, respectively, in good agreement with the experimental values 
of 3.467-i^ and 0.74 A i^. 

III. RESULTS AND DISCUSSION 

A. Hydrogen dissociation on the clean 7-U(100) surface 

The geometry and electronic properties of the clean 7-U(100) surface is firstly investi- 
gated. After geometry optimization, we find that the two topmost U layers relax significantly 
from the bulk values. The first-second interlayer is contracted by 26.4% and the second- 
third interlayer is expanded by 15.6%. The huge surface relaxation reflects that the surface 
electronic structure should be quite different from that in the bulk. The projected density of 
states (PDOS) around the Fermi energies are then calculated for bulk U and the 7-U(100) 
surface, and shown in Fig. 1(a). We can see that in the PDOS for bulk U, the occupied 
electronic states are itinerant without any localized peaks, which is in good agreement with 
previous theoretical studiesi^i^i. For the U atom in the 7-U (100) surface, the narrow peaks 
in the unoccupied / states disappear. This result supports the experimental observations 
in X-ray and ultraviolet photoelectron spectroscopy, and auger electron spectroscopy that 
localization effects are weak in U films^^. The PDOS for an U atom in a-U and in the 
a-U(OOl) surface is also shown in Fig. 1(b). One can see that there are more / electrons for 
the U atom in the a-U(OOl) surface than in a-U, which corresponds to the charge transfer 
from bulk to the surface atomic layer. Comparatively, such kind of charge redistribution is 
not seen for the 7-U(100) surface. 

After studying the geometry and electronic property of the 7-U(100) surface, we build our 
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models to calculate the two-dimensional (2D) potential energy surface (PES) cuts for H2 on 
the relaxed uranium surface. As shown in Fig. 2(a), there are three different high-symmetry 
sites on the clean U surface, respectively, the top, bridge (bri), and hollow (hoi) sites. At 
the bridge site, an adsorbed H2 has three different high-symmetry orientations, respectively 
along the x (i.e., [001]), y (i.e., [010]), and z (i.e., [100]) directions. Nevertheless, at the 
top and hollow sites, the x and y directions are degenerate, and the three high-symmetry 
orientations become the x, (i(i.e., [Oil]), and z directions. Herein, we use top-x, d, 2;, bri- 
X, 2, and hol-x, d, z to represent the nine high-symmetry channels for the adsorption of 
H2 on the 7-U(100) surface, respectively. We have also constructed several low-symmetry 
initial structures by rotating the H2 molecule in the XY ^ FZ, and XZ planes with small 
angles. Through PES calculations, we find that the energy barrier for H2 dissociation along 
the low-symmetry channels is always larger than along the high-symmetry ones. Similar 
results have also been obtained for the 02/Pb(lll) system where O2 adsorption also prefers 
the high-symmetry channels^^. Therefore, we will only discuss the obtained PES cuts along 
the high-symmetry channels. 

From our PES calculations, we find that there are no molecular adsorption states for H2 
on the 7-U(100) surface. The calculated 2D PES cuts along the top-x, hol-d, bri-y and bri-z 
channels are respectively listed in Figs. 3(a)-3(d). The PES cuts along the top- and hol-^ 
channels have similar energy distributions with the bri-z channel, and thus are not listed. 
Along the other bri-x, top-x, d and hol-x, d adsorption channels, the H2 molecule dissociates 
after overcoming small energy barriers. Only in the PES cut along the hol-c? channel, we 
find a local energy minimum after H2 dissociation, which is shown in Fig. 3(b). Along all 
the other dissociation channels, the H-H bond length is larger than 2.4 A in the atomic 
adsorption states after dissociation. So we do not see the local minima in the calculated 
PES cuts. The local minimum point in Fig. 3(b) corresponds to the adsorption state of two 
hydrogen atoms in the same surface uranium square hollow. After geometry optimization 
from this point, we find that the surface uranium atoms are distorted to lower down the 
free electronic energy. And in the stable adsorption state, the free energy of the adsorption 
system is 2.59 eV lower than that of an isolated H2 molecule plus a clean 7-U(100) surface. 

The most energetically favorable dissociation channel for H2 is along the bri-y channel on 
the 7-U(100) surface, which is found to be with no energy barriers. As shown in Fig. 3(c), 
the free electronic energy of the adsorption system drops by about 0.4 eV at the molecular 
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height of /iH2=0.90 A, when the hydrogen bond length c/h-h becomes 1.90 A. Our result 
that H2 molecules dissociate without any energy barriers on the 7-U(100) surface indicates 
different surface chemical properties of 7-U with a-U, because hydrogen dissociation on the 
a-U surface needs to overcome a small energy barrier of 0.08 eV^^. And as we will see in the 
following, the electronic interactions during H2 dissociation on the 7-U surface is different 
from that on the a-U surface. 

The PDOS evolution of H and U atoms along the bri-y channel is then analyzed to study 
the electronic interactions during the barrierless dissociation process of H2 on the 7-U(100) 
surface. We have chosen four points along the minimum energy dissociation path to calculate 
the PDOS, which are (/iH2=1.17 A, c/h-h=1-28 A), (/iH2=1.02 A, c;h-h=1-58 A), {h^^=0.90 
A, c?H-H=l-90 A), and (/ih2=0.84 A, dH-H=2.38 A) respectively. The obtained PDOS of the 
adsorption system at these points are listed in Figs. 4(a)-(d) respectively. As shown in Fig. 
4(a), at the molecular height of 1.17 A, a few electrons transfer from the molecular orbital 
of H2 to the unoccupied d states of U, which will be called as the charge donation process 
at following discussions. At the same time, some electrons also transfer back from U to H, 
forming a new peak in the hydrogen s states near 2.25 eV below the Fermi energy, and we 
will call it as a charge back-donation process. Since the change in free electronic energy is 
still negligible at the molecular height of 1.17 A , as shown in Fig. 3(c), the above electronic 
hybridizations are not strong. 

When the H2 molecule further approaches the 7-U surface, the electronic hybridizations 
become stronger, and the charge donation and back-donation become more obvious. As 
shown in Figs. 4(b)-4(d), more and more electrons transfer from s states of H to the 
unoccupied d states of U, and more and more electrons transfer from d and / states of U 
back to hydrogen. We can also see from Figs. 4(b)-4(d) that the electronic hybridization 
between H s states and U d states is always stronger than that with U / states. And the 
charge donation process only happens between H and d electronic states of U. Therefore, d 
electronic states of U play very important roles during the dissociation of H2 molecules on 
the 7-U surface. These results are quite different from the dissociation of H2 molecules on 
the a-U surface, where d electronic states of U is negligible, and the electronic hybridization 
happens between hydrogen s and uranium / electronic states^. 

To investigate the charge transfer for the dissociative adsorption of H2 on the 7-U(100) 
surface, we then calculate the atomic charges for the final state of the minimum energy 
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dissociation path, using the Bader topological method^. It is found that after dissociation, 
the two hydrogen atoms together gain 0.90 electrons, indicating that the ionic part of the 
H-U bonding plays a significant role during the dissociation process^^. Different from the 
dissociation of H2 on the a-U(OOl) surface, where the charge transfer happens only between 
the topmost U atoms and hydrogen atoms^^, here we find that the topmost and second 
atomic layer of the 7-U(100) surface both lose electrons to hydrogen atoms. In comparison 
with a bare relaxed surface, the topmost and second layer loses 0.58 and 0.34 electrons 
respectively. This result indicates that the electronic interaction between II2 and the 7- 
U(IOO) surface is not so localized as that between H2 and the a-U(lOO) surface, and also 
prove the participation of more itinerant d electronic states in interactions with hydrogen 
electrons. 

B. Hydrogen dissociation on the Nb-doped 7-U(100) surface 

Since the 7 phase of U is more stable after doping with 3d transition metal atoms. 
We here also investigate the influences of Nb-doping on the dissociation of H2 molecules 
on the 7-U(100) surface. Previous studies have already revealed that doped Nb atoms 
thermodynamically prefer to substitute U atoms in 7-U, rather than occupy octahedral or 
tetrahedral vacancies^- . Thus we only consider the substitutional doping of Nb atoms on 
the 7-U (100) surface. 

Firstly, we do geometry optimization for the Nb-doped U surfaces, with the Nb atom 
in the topmost, second, and third layer respectively. To simplify our discussions, we will 
call them as the UNbl, UNb2, and UNb3 surfaces in the following. The relaxed surface 
structures of the UNbl and UNb2 surfaces from the top view are shown in Figs. 2(b) and 
2(c) respectively. The adsorption sites of II2 molecules on them are the same as that on 
the clean 7-U(100) surface as depicted in Fig. 2(a). After geometry optimizations, the 
uranium atoms in the doping layer are no longer in the same plane with the doped Nb 
atoms. For example for the UNbl surface, the z coordinate of the Nb atom is 0.16 A larger 
than its nearest U atoms, and 0.09 A smaller than its next neighboring U atoms. The 
relative surface relaxation of Nb-doped 7-U surfaces can be calculated by averaging the z 
coordinates of the four atoms in the same layer. In this way, the relative relaxation (i.e., 
Adij/do with dij and c/q to be the interval between the ith and jth atomic layer, and the 
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lattice interval in bulk 7-U) is calculated and summarized in Table I for the UNbl, UNb2, 
and UNb3 surfaces. Prom the results hsted in Table I, we can see that the UNbl surface 
has much smaller surface relaxations than the clean 7-U (100) surface, indicating that their 
surface electronic structures are different. As the Nb atom is doped deeper, the surface 
relaxations tend to approach to the values of undoped 7-U (100) surface. 

The PDOS for the UNbl, UNb2, and UNb3 surfaces are also calculated and shown in 
Figs. 5(b)-5(d), in comparison with the PDOS of the undoped 7-U(100) surface shown in 
Fig. 5(a). We can see that the electronic states around the Fermi energies are contributed 
by both the d electrons of Nb, and d, / electrons of U. Therefore, one can expect that the 
surface interactions with atoms or molecules should be influenced by the doped Nb atoms, 
especially for the UNbl surface where the Nb atom is at the topmost layer. The PES for 
H2 on the UNbl surface is then calculated to study the influences. 

The calculated 2D PES cuts for H2 molecules on the UNbl surface along the top-x, 
hol-d, bri-y and hri-z channels are hsted in Figs. 6(a)-6(d) respectively. One can see great 
similarities in the energy distributions with that on the clean 7-U (100) surface. Firstly, along 
the bri-, top-, and ho\-z channels, H2 molecules can hardly adsorb or dissociate. Secondly, the 
most energetically favorable dissociation path is along the hh-y channel, which is without 
any energy barriers. And at last, the H2 dissociation along the other channels needs to 
overcome small energy barriers. Nevertheless, because of the introduction of a surface Nb 
atom, there are also some new features. As shown in Figs. 3(a) and 6(a), the dissociation 
energy barrier is larger on top of the Nb atom than on top of a U atom of the undoped 
U surface. And for the most energetically favorable dissociation path, we see that the free 
electronic energy lowers down much more quickly next to the Nb atom than next to a U 
atom of the undoped U surface, as shown in Figs. 3(c) and 6(c). These results indicate that 
the surface electron distribution changes after doping with a Nb atom. The larger energy 
reduction for H2 dissociation on the UNbl surface also suggests that surface doping of Nb 
atoms reinforces H2 dissociation, instead of hindering it. 

The electronic interactions between H2 and the UNbl surface are then studied for the 
most energetically favorable dissociation channel bri-y. The electronic evolution is analyzed 
by calculating the PDOS of H and U, Nb atoms at such four points: (/iH2 = l-50 A, (iH-H=0.84 
A), (/iH2=0.90 A, (iH-H=1.30 A), (/iH2=0.80 A, (iH-H=1.53 A), and {hn^^0.70 A, ciH-H=2.12 
A) along the dissociation path, which are shown in Figs. 7(a)-7(d) respectively. We can 
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see from Fig. 7(a) that at the molecular height of 1.50 A, no obvious electronic interactions 
happen between H and U, Nb atoms. When the molecular height lowers to be 0.90 A, some 
electrons transfer from H2 to the unoccupied Nb d electronic states, corresponding to the 
charge donation process, and some other electrons transfer back from d states of U and Nb 
to the antibonding orbital of H2, corresponding to the charge back-donation process. At 
the molecular height of 0.80 A, the charge donation and back-donation processes become 
more obvious, as shown in Figs. 7(e) and 7(f). We can also see from Fig. 7(e) that the 
charge donation mainly goes to the d orbitals around the Nb atom. In comparison with 
the electronic interaction between H2 and the undoped U surface, one can see that the d 
electronic states of Nb largely participate in both the charge donation and back-donation 
processes. And the stronger electronic hybridizations between hydrogen s and Nb d states 
lead to the fact that the free electronic energy lowers down more quickly on the doped 7-U 
surface. 

IV. CONCLUSIONS 

In conclusion, we have systematically studied the dissociation of H2 molecules on the 
clean and Nb-doped 7-U(100) surfaces. We find that both of the two reactions are with no 
energy barriers. On the clean 7-U (100) surface, there are electronic interactions between 
hydrogen electrons and U d electrons, which cause that H2 molecules dissociate without 
any energy barriers along the bri-y channel. This mechanism is quite different from the 
H2 dissociation on the a-\J surfaces, where / electrons of U, instead of d electrons take 
part in the electronic interactions. After surface doping of a Nb atom, we find that not 
only the H2 dissociation is with no energy barriers, but also the dissociated hydrogen atoms 
bond stronger with the surface. The d electronic states of Nb participate in the electronic 
interactions with hydrogen, and causes the larger free energy reduction. 
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TABLE I: The relative surface relaxation for the clean 7-U(100) surface, and the UNbl, UNb2, 
UNb3 surfaces, dij represents for the interval between the ith and jth atomic layer of each surface, 
and do is the lattice interval along the [100] direction of bulk 7-U. 



surfaces 


7-U(100) 


UNbl 


UNb2 


UNbS 


A(i23/c^0 


26.4% 


17.4% 


23.8% 


29.9% 


Adu/do 


15.6% 


7.0% 


21.9% 


16.0% 


Ad45/do 


9.8% 


0.7% 


10.4% 


11.1% 
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Fig.l (Color online) The projected density of states for a uranium atom in bulk and in 

the (100) surface of 7- (a) and a-U (b). The Fermi energies are set to be zero. 

Fig. 2 (Color online) (a) Top view of the clean 7-U(100) surface with the high-symmetry 
adsorption sites and high-symmetry H2 orientations depicted, (b) and (c) Top view of the 
UNbl [U(IOO) surface with a doping Nb atom in the topmost layer] and UNb2 [U(IOO) 
surface with a doping Nb atom in the second layer] surfaces, (d) Side view of the adsorption 
model for H2 on the clean 7-U(100) surface. Blue, grey, and pink balls represent for 
hydrogen, uranium, and niobium atoms respectively. The adopted supercell is depicted by 
dashed hues. 

Fig. 3 (Color online) The 2D PES cuts for the adsorption of hydrogen molecules along the 

(a) top-x, (b) hol-d, (c) bri-y, and (d) hri-z channels on the clean 7-U(100) surface. The to- 
tal energy of a free H2 molecule plus that of the 7-U(100) surface is set to be the energy zero. 

Fig. 4 (Color online) The projected density of states for the H2/7-U(100) adsorption 
system along the energetically most favorable dissociation path with the height of H2 center 
of mass to be 1.17 A (a), 1.02 A (b), 0.90 A (c), and 0.84 A (d). The Fermi energies are 
all set to be zero, (e) and (f) The partial charge density distributions for the two energy 
peaks denoted in (c). Blue and grey balls represent for hydrogen and uranium atoms 
respectively. The yellow dots represents for the isosurface of partial charge density. 

Fig. 5 (Color online) (a) The projected density of states for the clean 7-U(100) surface. 

(b) , (c), and (d) The projected density of states for the Nb doped 7-U(100) surface, with 
the Nb atom at the topmost, second, and third layer. The Fermi energies are all set to be 
zero. 

Fig. 6 (Color online) The 2D PES cuts for the adsorption of hydrogen molecules along the 
(a) top-x, (b) hol-d, (c) hri-y, and (d) bri-2; channels on the UNbl surface. The total energy 
of a free H2 molecule plus that of the Nb doped 7-U(100) surface is set to be the energy zero. 
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Fig. 7 (Color online) The projected density of states for the adsorption system of H2 on 
the Nb-doped 7-U(100) surface, along the energetically most favorable dissociation path 
with the height of H2 center of mass to be 1.50 A (a), 0.90 A (b), 0.80 A (c), and 0.70 
A (d). The Fermi energies are all set to be zero, (e) and (f) The partial charge density 
distributions for the two energy peaks denoted in (c). Blue, grey, and pink balls represent 
for hydrogen, uranium, and niobium atoms respectively. The yellow dots represents for the 
isosurface of partial charge density. 
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